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Abstract. - Recent studies on the Portevin - Le Chatelier effect report an intriguing crossover 
phenomenon from a low dimensional chaotic to an infinite dimensional scale invariant power 
law regime in experiments on CuAl single crystals and AlMg polycrystals, as a function of 
strain rate. We devise a fully dynamical model which reproduces these results. At low and 
medium strain rates, the model is chaotic with the structure of the attractor resembling the 
reconstructed experimental attractor. At high strain rates, power law statistics for the magni- 
tudes and durations of the stress drops emerge as in experiments and concomitantly, the largest 
Lyapunov exponent is zero. 



The Portevin-Le Chatelier (PLC) effect, discovered at the turn of the last century H, 
is a striking example where collective behaviour of defects leads to complex spatio-temporal 
patterns The PLC effect manifests itself as a series of serrations on the stress-strain 

curves when samples of dilute alloys are deformed under constant strain rate, e a (actually con- 
stant pulling speed). The effect is observed only in a window of strain rates and temperatures. 
Each stress drop is associated with the formation and often the propagation of a dislocation 
band. In polycrystals, at low e a , the randomly nucleated type C bands with large stress drop 
amplitudes are seen. At intermediate strain rates, one finds the spatially correlated 'hopping' 
type B bands moving in a relay race manner with smaller stress drop amplitudes. At high 
strain rates, propagating type A bands with small amplitudes are observed. (In single crystals 
such a clear classification does not exist.) These different types of PLC bands are believed to 
represent distinct correlated states of dislocations in the bands. It is this rich spatio-temporal 
kj^j | dynamics that has recently attracted the attention of physicists as well ^® ■ Indeed, the PLC 
effect is a good example of slow-fast dynamics commonly found in many stick-slip systems 
such as frictional sliding j5) , fault dynamics || and peeling of an adhesive tape |7j . 

Recent efforts have shown that surprisingly large body of information about the nature of 



dynamical correlations is hidden in the stress-strain curves 10, More recently, an 

intriguing crossover phenomenon from a chaotic regime occurring at low and medium strain 
rates to a power law regime at high strain rates has been detected in experiments on the PLC 
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effect in Cu-10% Al singie crystals [j0| and Al-Mg polycrystals |l3|]. This suggests that the 
crossover itself is insensitive to the micro structure. The purpose of this paper is to extend a 
model for the PLC effect introduced earlier |l4|] to explain this crossover phenomenon. 

This crossover phenomenon is of interest in the larger context of dynamical systems as this 
is a rare example of a transition between two dynamically distinct states. Chaotic systems are 
characterised by the self similarity of the strange attractor and sensitivity to initial conditions 
quantified by fractal dimension and the existence of a positive Lyapunov exponent, respec- 
tively. In contrast, a class of spatially extended dissipative systems often exhibit a tendency 
to evolve to a marginally stable state, characterised by power law statistics for the events, 
under the action of slow external drive with no parameter tuning. (It must be emphasized 
that power laws can arise due to other mechanisms also Such a state is termed as self- 

organised criticality (SOC) by Bak et al Unlike the low dimensional nature of the chaotic 
attractor, SOC state is an infinite dimensional state. Large number of physical systems are 
known to exhibit SOC type features 

There are numerous models and experiments where either of these dynamical regimes have 
been detected 

To the best of our knowledge ref. |]12| is one of the two instances 
known where both these states are observed in one and the same system. The other example 
is in hydrodynamics where chaos is observed at low Rayleigh number and power law scaling 
regime, known as hard turbulence, is seen at high values f2(i|j . Moreover, both the hard 
turbulence in convection and the power law regime in the PLC effect are observed at high 
drive parameter values in contrast to most SOC systems jl7],[l^]. Thus, power laws in the 
PLC effect are closer to the turbulent regime in convection than other SOC systems. 

The microscopic origin of the PLC effect is due to the interaction of mobile dislocations 
with diffusing solute atoms and is referred to as dynamic strain aging (DSA), first suggested 
by Cottrell Q and later improved by others. (See for instance ref. ||, ||). At low strain 
rates (or high temperatures) the average velocity of dislocations is low, there is sufficient 
time for solute atoms to diffuse to dislocations and pin them (usually called as aging). Thus, 
longer the dislocations are arrested, larger will be the stress required to unpin them. When 
these dislocations are unpinned, they move at large speeds till they are pinned again. At 
high strain rates (or low temperatures), the time available for solute atoms to diffuse to the 
dislocations decreases and hence the stress required to unpin them decreases. Thus, in a range 
of strain rates and temperatures where these two time scales are typically of the same order of 
magnitude, the PLC instability manifests. The competition between the slow rate of pinning 
and sudden unpinning of the dislocations, at the macroscopic level translates into a negative 
strain rate sensitivity (SRS) of the flow stress as a function of strain rate which is the basic 
instability mechanism used in most phenomenological models 

The well separated time scales implied in the DSA is mimicked by the fast mobile, immobile 
and the intermediate 'decorated' Cottrell type dislocations in the dynamical model due to 
Ananthakrishna and coworkers [[u} . The basic idea of the model is that qualitative features 
of the PLC effect emerge from the nonlinear interaction of these few dislocation populations, 
assumed to represent the collective degrees of freedom of the system. The rate equations for 
dislocation densities are constructed based on known dislocation mechanisms. In spite of the 
idealised nature of the model, it is successful in explaining several generic features of the PLC 



effect, such as the emergence of the negative SRS |14,|2l|], the existence of critical strain for the 



onset of the PLC instability, and the existence of a window of strain rates and temperatures 
for the occurrence of the PLC effect (see [Q). One prediction of the model is that there is a 
range of strain rates where the PLC effect is chaotic p^] , subsequently verified by analysing 
experimental signals (l0[ |ll|, [l^, [l3| . The model has been studied in detail by our group and 
others including an extension to the case of fatigue p3|,p^,p5|. 
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Since the model is fully dynamical and it predicts chaos at intermediate strain rates found 
in experiments, it has the right ingredients for studying this crossover. Here, we outline the 
model in terms of scaled variables in the notation of Ref. p6[ . In addition, we introduce a 
spatial coupling arising out of the cross-slip mechanism, used earlier by others [^|,|^]. The 
model consists of densities of mobile, immobile, and Cottrell's type dislocations denoted by 
p m (x,t), pim(x,t) and p c (x,t) respectively, in the scaled form. The evolution equations are: 

d Pm , 2 , , , m D d 2 (^ ff (x) Pm ) 

= -b p m - pmPim + pim ~ apm + <P e ffPm + ~ 7^ > U) 



dpi 



bo(b Pm - PmPim ~ Pirn + ape), (2) 
C{p m ~ Pc)- (3) 



dt 

dp c 
dt 

The first term in eqn.(l) refers to the formation of locks and consequent immobilisation of 
two mobile dislocations, the second term to the annihilation of a mobile dislocation with 
an immobile one, and the third term to the remobilisation of the immobile dislocation due 
to stress or thermal activation. The fourth term represents the immobilisation of mobile 
dislocations due to solute atoms. Once a mobile dislocation starts acquiring solute atoms 
we regard it as the Cottrell's type dislocation p c . As they progressively acquire more solute 
atoms, they eventually stop, then they are considered as immobile dislocations pi m (loss term 
in eqn. (3) and the gain term in eqn. (2)). Alternately, the aggregation of solute atoms 
can be regarded as the definition of p c , ie., p c = f_ dt' p m (t')K(t — t'). For the sake of 
simplicity, we use a single time scale with K(t) = e~ ct . The convoluted nature of the integral 
physically implies that the mobile dislocations to which solute atoms aggregate earlier will 
be aged more than those which acquire solute atoms later ( see ref. |2l|]). The fifth term 
represents the rate of multiplication of dislocations due to cross-slip. This depends on the 

1 /2 

velocity of the mobile dislocations taken to be V m {4>) = ^Tff^ where 4>ef / = (0 — hp^ ) is the 
scaled effective stress, <j> the scaled stress, m the velocity exponent and h a work hardening 
parameter. The last term is a spatial coupling term arising out of the nonlocal nature of 
the cross-slip as argued below. Actually, the nature of the spatial coupling in the PLC 
effect has been a matter of debate [| . Within the scope of the model, a natural source 
of spatial coupling is the nonlocal nature of cross-slip as dislocations generated at a point 
spread over to the neighbouring elements. (Compatibility stresses between the slipped and 
the unslipped regions and long range interactions are other possible sources of coupling ||[||].) 
Let Ax be an elementary length. Then, the flux $(a;) flowing from x ± An; and out of x 
is given by $(x) + § [$(x + Ax) - 2§(x) + $(x - Ax)] where $(x) = p m (x)V m (x). Here p 
is the probability of cross-slip spreading into neighbouring elements. Expanding $(x ± Ace) 

and keeping the leading terms, we get p m V m + § 9 %™^ m ^ (Az) 2 . Further, cross-slip spreads 
only into regions of minimum back stress. Noting that the back stress is usually taken to 
result from the immobile dislocation density ahead of it, we use Ax 2 —< Ax 2 >= f 2 p~^. 
Here < . . . > refers to the ensemble average and f 2 is an elementary (dimensionless) length. 
Finally, a, bo and c are the scaled rate constants referring, respectively, to the concentration 
of solute atoms slowing down the mobile dislocations, the thermal and athermal reactivation 
of immobile dislocations, and the rate at which solute atoms are gathering around the mobile 
dislocations. We note that the order of magnitudes of the constants have been identified in 
ref. Jl4,23|. These equations are coupled to the machine equation 

^M=d[e-j f p m (x,t)<p™ ff (x,t)dx], (4) 
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where e is the scaled applied strain rate, d the scaled effective modulus of the machine and the 
sample, and I the dimensionless length of the sample. (We reserve e a for the unsealed strain 
rate.) 

The PLC state is reached through a Hopf bifurcation. The domain of instability in e is 
10 < e < 2000, and that in other parameters is the same as in the original model, beyond which 
uniform steady states exist . Here, we use a = 0.8, bo = 0.0005, c = 0.08, d = 0.00006, m = 
3.0, h — 0.07 and D = 0.5. (The parameter values used here are essentially the same used in 
several of our earlier calculations. ^(|.) But the results discussed below hold true for a wide 
range of other parameters in instability domain including a range of values of D. 

These equations are discretised into M equal parts of width Al and solved numerically for 
Pm(jit), Pim(jjt), PcU't)? J = 1: 2, M , and </>(i). Due to the widely differing time scales, 
appropriate care has been exercised in the numerical solutions of these equations by using 
a variable step fourth order Runge-Kutta scheme with an accuracy of 10 -4 for all the four 
variables. The spatial derivative of p m is approximated by its central difference. The initial 
values of dislocation densities are taken to be their steady state values with a Gaussian spread 
along the length of the sample with a Gaussian spread along the length of the sample, since the 
long term evolution of the system is essentially independent of the initial values. Now consider 
the boundary conditions. Since the sample is strained at the grips, we choose Pi m (j, t) at j = 1 
and N to be three orders of magnitude more than rest of the sample and p m (j, t) — p c {j, t) =0 
at j = 1 and N, as p m and p c cannot evolve into the grips . 

We first recall the relevant experimental results on the crossover phenomenon and compare 
them with those from the model. Plots of two experimental stress-strain curves corresponding 
to the chaotic and SOC regimes of applied strain rates are shown in Fig. 1. In ref. jL2j, the 
chaotic nature of the stress-strain curve shown in Fig. la, was demonstrated by showing the 
existence of a finite correlation dimension using the standard Grassberger-Procaccia algorithm 
, and the existence of a positive Lyapunov exponent Both these methods involve 

embedding the time series in a higher dimensional space using time-delay technique [ p"9| . (In 
addition, surrogate data analysis was also carried out in Jl2[.) The correlation dimension, v, 
of the experimental attractor was found to be 2.3. Then, the number of degrees of freedom 
required for the description of the dynamics of the system is given by the minimum integer 
larger than v + 1 which is four in this case, consistent with that used in the original model. 
The geometrical interpretation of these degrees of freedom is that it is the subspace to which 
the trajectories are confined. This dimension can also be obtained by an alternate method, 
called the singular value decomposition psj , which has an additional advantage of allowing 
the visualisation of the strange attractor. This method has been applied to the PLC time 
series earlier In this method, the trajectory matrix is constructed and the eigen values 
of the covariance matrix are calculated. For the time series in Fig. la, we find that the 
relative strength of the fourth eigen value drops more than two orders of magnitude compared 
to the first and changes very little beyond the fourth eigen value. Thus, we estimate the 
dimension of the experimental attractor to be four, consistent with that obtained from the 
correlation dimension. Using the first three principal directions of the subspace Q; i = 1 to 3, 
we have reconstructed the experimental attractor in the space of specifically chosen directions 
C\ — C2,Cs and C\ to permit comparison with the model. This is shown in Fig. 2a for the 
experimental time series at e a = 1.7 x 10~ 5 s _1 . This can be compared with the strange 
attractor obtained from the model in the space of pmipim an d Pc (at an arbitrary spatial 
location, here j = 50 and N — 100) shown in Fig. 2b for e = 120 corresponding to the mid 
chaotic region (see below). Note the similarity with the experimental attractor particularly 
about the linear portion in the phase space (Fig. 2a) identified with the loading direction in 
Fig. la. Note that the identificaton of the loading direction is consistent with the absence 
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of growth of p m . Since the existence of a positive Lyapunov exponent is a confirmatory test 
of chaos, we have calculated both the largest Lyapunov exponent (LLE) and the spectrum of 
Lyapunov exponents ]l9| directly from the model. Since the magnitude of the LLE and that 
of the largest Lyapunov exponent calculated from the spectrum agree, we have shown the 
average LLE (obtained over 15000 time steps after stabilisation) in Fig. 3a as a function of 
the strain rate for N = 100. The LLE becomes positive around e = 35 reaching a maximum 
at e = 120, and practically vanishing around 250. (Periodic states are observed in the interval 
10 < e < 35.) For e > 250, the dispersion in the value of the LLE is ~ 5 x 10 -4 which is the 
same order as the mean. Thus, the LLE is taken to vanish beyond e = 250 as is expected of 
the power law regime p9[ |. 

At high strain rates beyond i ~ 280, the stress-time series obtained from the model shows 
no inherent scale in the magnitudes of the stress drops as in the case of experimental time series 



at high strain rates |12,13j. We have analysed the distributions for the stress drop magnitudes 
Atfi and their durations At. The distribution of stress drop magnitudes, D(A<fi), shows a power 
law D(A(f>) ~ A<j)~ a . This is shown in Fig. 3b(o) along with the experimental points (•) 
corresponding to e a = 8.3 x 10 s" 1 . Clearly both experimental and theoretical points show 
a scaling behaviour with an exponent value a m 1.1. The distribution of the durations of the 
drops D(At) ~ At" 13 also shows a power law with an exponent value (5 « 1.3. The conditional 
average of A<fi denoted by < Acj) > c for a given value of At behaves as < A(f> > c ~ At 1 ^ with 
x w 0.65. The exponent values satisfy the scaling relation a = x(f3 — 1) + 1 quite well. Since 
the basic cause of the stress-drops is the growth of mobile dislocation density during this 
period, one can look at the scaling behaviour of the total density in the sample at a given 
time using p m (t) = J p m (x,t)dx. Let Ap m (t) denote the increase in p m (t) occurring during 
the intervals of the stress drops. Then, one should expect that the statistics of Ap m also to 
exhibit a power law, ie., D(Ap m ) ~ Ap~ 7 , with the same exponent value as a. A plot of 
D(Ap m ) for N — 300 is shown in Fig. 3b ({>). The extent of the power law regime is nearly 
two orders with 7 » 1.1, same as a. The large bump at high values is due to the effect of finite 
size of the system as in many models || and in particular, here it is due to high levels of stress 
at the grips. Noting that dislocation bands cannot propagate into regions of high stresses, it 
is clear that the edges cause distortions in the otherwise smoothly propagating bands leading 
to large changes in p m (t). Typically, the influence of the edges are felt by the band when it 
is 20 sites away. Increasing N from 100 to 300, increases the scaling regime by half a decade 
and the peak of the bump reduces from 700 to 500, thus indicating the influence of the finite 
size of the system. 

To sum up, the present model exhibits chaotic dynamics at low and medium strain rates, 
and SOC dynamics at high e. These distinct dynamical states and the crossover clearly emerge 
due to the inclusion of nonlinearities in the form of basic dislocation mechanism and a spatial 
coupling. As for the nature of dynamics in different regions of strain rate, one can perhaps 
anticipate the emergence of the spatio-temporal chaotic regime |30| , as the original model 
exhibits chaos. However, the emergence of SOC dynamics needs some explanation. Recently, 
we have shown that the upper limit of the PLC effect in the original model is a result of a 
reverse Hopf bifurcation at high strain rates [^6|. This implies that the average amplitude 
of the stress drops decreases as a function of strain rate as in experiments. This means 
that the average stress level (in time) is roughly constant with small fluctuations around the 
mean implying that e essentially balances the total plastic strain rate. The picture essentially 
remains unaltered when spatial coupling is introduced. Thus, in this regime, the state of the 
system is critically poised as in any SOC system. Recently, we have shown that the geometry 
of the slow manifold has a bent structure and the 'fold line' on the slow manifold corresponds 
to the threshold value of stress for unpinning the dislocations. (See figs. 4 and 5 in ref. |p3 
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and also for further details.) Interestingly, for the SOC regime, we find that most of the spatial 
elements are literally on the 'fold line' corresponding to the marginally stable state. 

The dynamical origin of the scaling regime in our model is similar to that dealt by Gil 
and Sornette (3^] using a subcritical Hopf bifurcation. More recently, a dynamical analysis 
of Zhang's model of SOC has been reported as well p^] . However, it must be pointed out 
that in the present case, the scaling occurs at high drive values in contrast to many SOC type 
models. In this sense, the dynamical regimes found in our model are similar to that observed 
in thermal convection of a box of helium gas, namely, periodic states — ► chaos — > power law 
hard turbulence regime [ p0| as a function of Rayleigh number. (Actually, chaos and hard 
turbulence are separated by soft turbulence.) Lastly, to the best of our knowledge, this is the 
first model which exhibits such a crossover not just in the context of the PLC effect, but as a 
general crossover between two distinct dynamical regimes. 

* * * 

We thank Prof. Neuhauser for supplying the experimental data. This work is supported 
by Department of Science and Technology, Grant No: SP/S2/K-26/98, New Delhi, India. 

REFERENCES 

[1] Le Chatelier F., Rev. de Metallurgie, 6 (1909) 914; Portevin A. and Le Chatelier F., 

C.R. Acad. Set. Pans, 176 (1923) 507. 
[2] Viewpoint Set No 21, Propagative Instabilities. Scripta. Metall., 29 (1993) 1147. 
[3] Kubin L.P., Fressengeas C. and Ananthakrishna G., Collective Behaviour of Dislocations, 

edited by Nabarro F.R.N. , Dislocations in Solids, Vol. 11 (North-Holland, Amsterdam) 2002. 
[4] Cottrell A.H., Dislocations and Plastic Flow in Crystals (University Press, Oxford) 1953. 
[5] Persson B.N.J, and Tosatti E., Physics of Sliding Friction (Kluwer Academic Publishers, 

Dordrecht) 1996. 

[6] Carlson J.M. and Langer J.S., Phys. Rev. Lett., 62 (1989) 2632; Phys. Rev. A., 40 (1989) 
6470. 

[7] Maugis D. and Barquins M., Adhesion, edited by Allen K.W., Vol. 12 (Elsevier, London) 
1988, p. 205. 

[8] D'Anna G. and Nori F., Phys. Rev. Lett, 85 (2000) 4096. 

[9] Franklin S.V., Mertens F. and Marder M., Phys. Rev. E., 62 (2000) 8195. 
[10] Ananthakrishna G., et al.„ Scripta. Metall, 32 (1995) 1731. 
[11] Noronha S.J., et al, Int. J. of Bifurcation and Chaos, 7 (1997) 2577. 
[12] Ananthakrishna G., et al, Phys. Rev.E., 60 (1999) 5455. 
[13] Bharathi M.S., et al, Phys. Rev. Lett, 87 (2001) 165508. 
[14] Ananthakrishna G. and Valsakumar M.C., J. Phys. D., 15 (1982) L171. 
[15] See for instance chapter 14 of Sornette D., Critical Phenomena of Natural Sciences: Chaos, 

Fractals, Self-Organization and Disorder: Concepts and Tools (Springer, Heidelberg) 2000. 
[16] Bak P., Tang C. and Wiesenfeld K., Phys. Rev. Lett, 59 (1987) 381 -Phys. Rev. A., 38 
(1988) 364. 

[17] Jensen H.J., Self- Organised Criticality (Cambridge University Press, Cambridge) 1998. 
[18] Bak P., How Nature Works (Springer - Verlag, NewYork) 1996. 

[19] Abarbanel H.D.I. , Analysis of Observed Chaotic Data (Springer- Verlag, New York) 1996. 

[20] Heslot F., Castaing B. and Libchaber A., Phys. Rev. A., 36 (1987) 5870. 

[21] Rajesh S. and Ananthakrishna G., Phys. Rev. E., 61 (2000) 3664. 

[22] Ananthakrishna G. and Valsakumar M.C., Phys. Lett., A95 (1983) 69. 

[23] Bekele M. and Ananthakrishna G.,Phys. Rev. E., 56 (1997) 6917. 

[24] Glazov M.V., et al, Appl. Phys. A, 64 (1997) 373. 



M.S. BHARATHI and G. ANANTHAKRISHNA: CHAOTIC AND POWER LAW STATES IN THE PORTEVIN-LE CHATELIER EFFECT7 



[25] Zaiser M., et al, Comp. Mater. Science, 5 (1999) 35. 

[26] Rajesh S. and Ananthakrishna G., Physica D, 140 (2000) 193. 

[27] P. Grassberger and I. Procaccia, P/w/sica D, 9 (1983) 189. 

[28] Broomhead D. and King G., Physica D, 20 (1987) 217. 

[29] Chen K., Bak P. and Obukhov S.P. ,Phys. Rev. A., 43 (1991) 625. 

[30] See Cross M.C. and Hohenberg P.C., Rev. Mod. Phys., 65 (1993) 851 for a discussion on 

the transition from chaos to spatio-temporal chaos. 
[31] Gil L. and Sornette D., Phys. Rev. Lett, 96 (1996) 3991. 

[32] Blanchard Ph., Cessac B. and Kruger T., J. Stat. Phys., 88 (1997) 307; J. Stat. Phys., 98 
(2000) 375; Cessac B., Blanchard Ph. and Kruger T., Phys. Rev. E., 64 (2001) 016133. 




Fig. 1 - Experimental stress-time series corresponding to (a) chaotic state at e a = 1.7 x 10 5 s 1 and 
(b) SOC state at e a = 8.3 x 10~ 5 s~\ 




Fig. 2 - (a) Reconstructed experimental attractor (b) Attractor from the model for N = 100, j = 50. 




Fig. 3 - (a) . The largest Lyapunov exponent of the model, (b) Distributions of the stress drops from 
the model (o), from experiments (•) and Ap m (0) from the model. Solid lines are guide to the eye. 



